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Abstract 


An  Image  Feature  Based  Space  Time  Processing  (IFSTP)  algorithm  is  introduced  to 
effectively  detect  moving  ground  targets  in  clutter  and  jamming  via  airborne  radar.  This 
new  approach  exploits  the  distinct  image  features  of  targets  and  interference  signals  in 
the  angle-Doppler  domain.  An  image  segmentation  algorithm,  referred  to  as  Region 
Growing  (RG),  extracts  targets  and  interference  features  in  the  angle-Doppler  domain, 
and  an  innovative  block-size  detection  algorithm  discriminates  between  moving  targets 
and  interference  based  on  the  extracted  image  features.  The  proposed  IFSTP  algorithm  is 
particularly  suitable  for  detecting  ground  moving  targets  in  highly  non-homogeneous 
clutter  environments,  without  any  requirement  for  clutter  covariance  estimation. 


PART  I:  IMAGE  FEATURE-BASED  SPACE  TIME 
PROCESSING  FOR  GROUND  MOVING  TARGET 
DETECTION 

1.  Introduction 

Space  Time  Adaptive  Processing  (STAP)  has  been  widely  used  in  advanced 
airborne  and  space-based  radar  systems  to  detect  ground  moving  targets  embedded  in 
heavy  clutter  and  jamming  environments  [1-5].  Joint  space-time  processing  is  necessary 
for  ground  clutter  rejection  because  the  ground  clutter  couples  between  the  spatial 
domain  and  the  time  domain  due  to  platform  motion.  However,  successful 
implementation  of  STAP  requires  accurate  real-time  estimation  of  the  clutter  covariance 
matrix.  The  “training”  data  used  for  clutter  estimation  is  normally  obtained  by  sampling 
secondary  cells  that  are  spatially  adjacent  to  the  primary  detection  cell,  with  the 
assumption  that  the  clutter  in  the  primary  and  secondary  cells  is  statistically  Independent 
and  Identically  Distributed  (IID)  [1-3].  Moreover,  the  accuracy  of  the  clutter  estimate 
depends  upon  the  number  of  available  IID  secondary  data  [1].  Since  ground  clutter  is 
non-homogeneous,  finding  sufficient  IID  secondary  data  for  detection  processing  poses 
the  most  serious  challenge  to  successfully  implementing  STAP  algorithms.  Even  though 
rank-reduction  approaches  such  as  partially  adaptive  STAP  or  subspace-based  STAP  can 
lessen  the  requirements  on  secondary  data  [3,  4],  target  detection  using  STAP  is  still 
contingent  upon  the  availability  of  sufficient  IID  secondary  data.  In  some  extreme 
situations,  no  IID  secondary  data  may  be  available,  which  severely  limits  the 
effectiveness  of  traditional  STAP  methods.  This  paper  develops  a  novel  space-time 
processing  algorithm,  termed  Image  Feature  Based  Space-Time  Processing  (IFSTP),  for 
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effective  ground  moving  target  detection  in  non-homogeneous  clutter  environment 
without  using  any  secondary  data. 
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2.  Image  Feature-Based  Space-Time  Processing  (IFSTP) 

Consider  an  airborne  sensor  system  that  collects  echo  data  in  the  space-time  domain 
by  transmitting  multiple  coherent  pulses  and  receiving  data  from  each  element  of  an 
antenna  array.  STAP  suppresses  interference  by  “whitening”  interference  signals  and 
further  integrates  target  signals  through  two-dimensional  matched  filtering  in  the  time- 
space  domain  [5],  The  proposed  IFSTP  technique,  however,  transfonns  the  collected 
space-time  data  directly  into  the  angle-Doppler  domain;  hence,  both  target  signal  and 
interference  are  coherently  integrated  through  the  transform.  The  discrimination 
processing  of  targets  and  interference  signals  for  detection  processing  is  performed  based 
on  their  different  features  in  the  angle-Doppler  plane  without  explicitly  filtering  out 
interference  signals. 

Figure  1  shows  a  typical  target-and-interference  image  of  a  detection  cell  in  the 
angle-Doppler  plane  for  an  airborne  radar  system. 


Angle 


Figure  1:  Typical  target  and  interference  scene  image  on  the  angle-Doppler  plane  for 

airborne  radar  systems 
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If  the  number  of  transmitted  pulses  during  a  Coherent  Processing  Interval  (CPI)  and 
the  number  of  array  elements  are  large  enough,  the  target  signal  becomes  a  concentrated 
point  in  the  transform  domain.  Jamming  is  strongly  correlated  in  the  angle  domain,  but 
completely  random  in  the  Doppler  domain.  Consequently,  as  shown  in  Fig.  1,  jammer 
signals  in  the  angle-Doppler  plane  are  straight  lines  parallel  to  the  Doppler  axis.  Due  to 
the  coupling  between  the  space  and  time  domains  through  antenna  sidelobes,  ground 
clutter  is  linearly  extended  in  both  angle  and  Doppler  domains.  As  a  result,  clutter 
normally  manifests  itself  as  a  tilted  ridge  or  even  multiple  parallel  ridges  depending  on 
the  ratio  of  the  platform  speed  to  the  radar  Pulse  Repetition  Frequency  (PRF).  Thennal 
noise,  not  shown  in  Fig.  1,  is  uncorrelated  in  both  the  space  and  time  domains,  and  thus 
statistically  uniformly  spread  in  the  angle-Doppler  plane  through  the  2-D  Fourier 
transform.  Therefore,  in  the  angle-Doppler  domain,  the  structure  of  target  signals  is 
conspicuously  different  from  that  of  ground  clutter  or  jamming,  i.e.  target  signals  are 
concentrated  while  interference  is  extended.  Note  that  moving  target  signals  generally  do 
not  overlap  with  clutter  on  the  angle-Doppler  plane  because  of  their  different  Doppler 
frequencies.  Hence,  based  on  the  above  observations,  we  develop  the  IFSTP  algorithm 
for  automatic  discrimination  of  moving  target  signals  from  interference  in  the  angle- 
Doppler  domain. 

The  first  step  of  IFSTP  is  to  transform  radar  data  collected  in  the  space-time 
domain  to  an  image  in  the  angle-Doppler  domain  through  a  two-dimensional  Fourier 
transformation.  Subsequently,  a  clamping  processing  is  applied  to  all  pixels  of  the 
transformed  images  to  remove  white  noise.  The  value  of  a  pixel  is  set  to  zero  if  its 
magnitude  is  smaller  than  a  pre-defined  threshold;  otherwise,  the  pixel  remains 
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unchanged  and  is  labeled  as  a  non-zero  pixel.  The  threshold  level  T  used  in  the  clamping 
processing  is  derived  from  the  standard  deviation  of  the  transfonned  white  noise, 

T  =  C-an  (1) 

where  C  is  a  constant  that  normally  is  chosen  to  be  between  2  and  3,  and  crn  is  the 

standard  deviation  of  the  white  noise,  which  is  mainly  generated  in  the  radar  receiver  [6]. 

Following  the  clamping  processing,  the  remaining  non-zero  pixels  are  either 
target  signals  or  interference  (clutter  or  jamming).  Feature  extraction  from  the  clamped 
image  is  carried  out  by  clustering  non-zero  pixels  into  pixel  blocks  consisting  of 
consecutively  connected  non-zero  pixels  through  an  image  segmentation  algorithm  called 
region  growing  [7,  8].  The  growing  of  a  pixel  block  is  implemented  by  comparing  the 
ranges  between  the  pixel  block  and  nearby  pixels  with  a  pre-defined  growing  distance  D. 
The  range  R  between  a  pixel  q  and  a  pixel  block  consisting  of  L  pixels  {pQ ,  px , A  pL  }  is 
defined  as  the  minimum  distance  between  pixel  q  and  any  pixel  in  the  block: 

R  =  min||/> . ,  g||,  1  <  i  <  L  (2) 

The  region  growing  procedure  we  developed  for  IFSTP  is  detailed  as  follows: 

Step  1:  Choose  any  non-zero  pixel  p'0  as  an  initial  element  of  block  f?„  i.e. 

Bi={p'0},i  =  0 

Step  2:  Grow  B,  by  merging  all  non-zero  pixels  p\,p\, A  , p'L  whose  ranges  to  Bj 
are  not  larger  than  the  pre-defined  growing  distance  D,  i.e. 

Bi  ^Bi  u{p[,p'2,A  ,plL)  (3) 

Step  3:  Repeat  Step  2  for  the  updated  Bt  until  no  more  non-zero  pixels  are  available 
for  merging  into  Bj. 
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Step  4:  Setting  /'<=/'  + 1  and  repeating  Steps  1-3  until  all  non-zero  pixels  are 
merged  into  blocks. 

Step  5:  The  final  result  consists  of  a  group  of  non-zero  pixel  blocks  B0,Bl,B2,L  . 

With  segmentation  processing,  the  image  becomes  a  collection  of  pixel  blocks  that  are 
either  target  or  interference.  Target  detection  for  a  pixel  block  is  based  on  pixel 
concentration  level,  which  may  be  measured  by  a  metric  called  block  size.  The  block 
size  Sk  of  a  pixel  block  Bk  is  defined  as  the  maximum  distance  of  any  two  pixels  inside 
the  block,  i.e. 

Sk  =  maxllpf ,  p].  II,  1  <  <  H  (4) 

ij  ii  ii 

where  H  is  the  number  of  pixels  in  the  block  Bk- 

Based  on  the  block  size  definition  in  (4),  the  following  detection  criterion  is 

employed  to  detennine  whether  a  block  Bk  with  size  Sk  in  the  angle-Doppler  domain  is  a 

target  signal  or  interference: 

SK  <  (3  =>  The  block  is  a  target 
SK  >  /3  =>  The  block  is  interfence 

where  the  detection  threshold  ft  is  detennined  from  the  image  resolution,  the  number  of 
pulses  during  a  CPI,  and  the  number  of  antenna  array  elements.  Because  an  ideal  point 
target  signal  could  occupy  up  to  four  resolution  cells,  the  minimum  value  of  fi  is  chosen 
to  be  2 r,  where  r  is  the  image  resolution.  With  further  consideration  of  such  block¬ 
expanding  factors  as  windowing,  limited  numbers  of  coherent  pulses,  and  antenna  array 
elements,  the  detection  threshold  /?is  selected  as: 

/?  =  (4  ~  8)  •  r  (6) 

The  detection  process  in  (5)  is  repeated  for  all  pixel  blocks  in  the  transfonned  image. 
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3.  IFSTP  Processing  Results 


We  applied  the  proposed  IFSTP  algorithm  to  radar  signals  generated  from  a 
simulated  airborne  radar  system.  The  radar  frequency  is  450MFIz,  the  platfonn  height  is 
9000m  and  the  range  is  130km.  A  linear  side-looking  antenna  array  with  32  elements 
and  X/2  element  spacing  is  employed.  There  are  32  coherent  pulses  in  a  CPI,  and  the 
values  of  SNR,  CNR  and  JNR  per  pulse  are  OdB,  40dB  and  60dB,  respectively.  The 
amplitude  of  the  received  target  signals  during  a  CPI  is  assumed  to  be  constant. 

Figure  2  shows  the  Fourier  transfonn  of  radar  echoes  containing  a  target,  jammer, 
and  clutter.  The  Doppler  frequency  of  the  target  is  -90FIz  and  its  azimuth  is  10°.  The 
jammer  is  located  at  an  angle  of  -30°.  The  radar  platfonn  velocity  is  50m/s.  As  a  result, 
the  slope  of  the  clutter  ridge  is  one  in  the  angle-Doppler  plane. 


-1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

SIN(PHI) 


Figure  2:  2-D  Fourier  transformed  image  on  the  Doppler-angle  plane  for  a  target  with 
Doppler  frequency  -90Hz  and  platform  velocity  50m/s. 
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Because  of  aliasing  effects  from  the  much  stronger  clutter  antijamming,  the  target  signal 
is  not  recognizable  in  Fig.  2.  If  an  85-dB  Chebyshev  window  is  applied  to  both  Doppler 
and  angle  domains,  and  the  transformed  image  is  subsequently  clamped  with  the 
threshold  given  in  (1),  the  target  signal,  jammer  and  clutter  ridge  become  clearly  visible, 
but  they  are  also  artificially  widened  due  to  the  windowing  processing.  Because  of  the 
clutter  ridge  expansion  and  7J2  array  element  spacing,  some  portions  of  the  grating  lobes 
of  the  clutter  ridge  appear  at  the  upper-left  and  lower-right  comers  of  the  image. 
However,  they  can  easily  be  eliminated  as  non-targets  because  they  are  located  at  the 
edges  of  the  image.  A  more  elegant  solution  is  to  allow  the  region  growing  algorithm  to 
form  pixel  blocks  by  warping  over  the  edges  of  the  image,  which  includes  the  false  target 
pixels  in  the  interference  blocks. 

For  extraction  of  target  and  interference  features,  the  minimum-distance  based 
region  growing  algorithm,  developed  in  Section  II,  is  applied  to  the  image  in  Fig.  3. 


-1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

SIN(PHI) 

Figure  3:  Clamped  Doppler-angle  image  with  85-dB  Chebyshev  window  tapering  for  radar 

echo  data  shown  in  Fig.  2. 
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The  minimum  distance  D  in  (2)  is  chosen  to  be  two  times  larger  than  the  image 
resolution.  Such  a  choice  is  essential  because  the  amplitude  of  interference  is  Rayleigh- 
distributed  and  there  are  some  null-points  in  the  sidelobe  of  the  antenna  array  pattern. 
The  non-zero  pixel  blocks  generated  using  the  region  growing  are  displayed  in  Fig.  4. 


(c)  (d) 

Figure  4:  Four  pixel  blocks  generated  from  the  image  in  Fig.  3  using  our  minimum- 
distance  based  region  growing  algorithm 


Finally,  the  block-size  based  target  detection  scheme  in  (5)  is  carried  out  for  all 
the  blocks  in  Fig.  4.  The  threshold  block  size  /?  used  for  target  detection  is  selected  as 
six.  The  detection  results  are  listed  in  Table  1. 
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Table  1:  Block  sizes  and  detection  results  for  the  generated  pixel  blocks  in  Fig.  4 


Block 

a 

b 

c 

d 

Block  Size 

43.8 

2.2 

1.4 

2.2 

Is  at  Image 
Edges? 

No 

Yes 

No 

Yes 

Target? 

No 

No 

Yes 

No 

Although  there  are  three  blocks  that  meet  the  target  size  criterion,  two  of  them  are 
at  the  image  edges  and  not  considered.  As  shown  in  Table  1,  block  c  is  the  only  target 
block  detected.  If  the  actual  target  is  located  at  the  image  edge  or  overlaps  with 
interference  signals,  a  slightly  different  radar  PRF  has  to  be  used  to  shift  the  moving 
target  signal  away  from  the  interference. 

Figure  5  is  a  generated  scene  that  includes  a  target  and  clutter,  but  no  jamming  in 
the  angle-Doppler  domain.  The  target  Doppler  frequency  is  -125  FIz  and  its  azimuth 
angle  is  0°. 


-1  -0.8  -0.6  -0.4  -0.2  0  0.2  0.4  0.6  0.8  1 

SIN(PHI) 


Figure  5:  2-D  Fourier  transformed  image  on  the  Doppler-angle  plane  for  a  target  with 
Doppler  frequency  -125  Hz  and  platform  velocity  100  m/s. 
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Figure  6  is  a  post-processing  result  of  the  image  in  Fig.  5,  when  a  Chebyshev  window 


and  clamping  are  applied  together. 
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Figure  6:  Clamped  Doppler-angle  image  with  Chebyshev  window  tapering  for  radar  echo 

data  shown  in  Fig.  5. 

With  the  velocity  of  the  radar  platform  increased  to  lOOm/s,  the  slope  of  the  clutter 
ridge  in  the  angle-Doppler  plane  becomes  two,  leading  to  multiple  parallel  clutter  ridges 
in  the  transformed  image.  Likewise,  using  the  region  growing  algorithm  the  target  and 
clutter  pixel  blocks  are  successfully  separated,  as  shown  in  Fig.  7. 
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Figure  7:  Four  pixel  blocks  generated  from  the  image  in  Fig.  6,  using  our  minimum- 
distance  based  region  growing  algorithm 


With  a  block-size  threshold  of  6,  the  block-size  based  detection  result  is  listed  in 
Table  2.  Block  c  is  correctly  detected  as  the  target. 

Table  2:  Block  sizes  and  detection  results  for  the  generated  pixel  blocks  in  Fig.  7 


Block 

a 

b 

c 

d 

Block  Size 

21.0 

36.4 

2.2 

20.1 

Is  at  Image 
Edges? 

No 

No 

No 

No 

Target? 

No 

No 

Yes 

No 
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Conclusions 


The  Image  Feature  Based  Space-Time  Processing  (IFSTP)  approach  is  developed  as 
an  alternative  approach  to  STAP  for  moving  ground  target  detection  in  clutter  and 
jamming.  Our  simulations  suggest  that  IFSTP  can  effectively  detect  moving  ground 
targets  based  on  distinguishing  features  of  targets  and  interference  in  the  angle-Doppler 
domain.  By  avoiding  the  requirement  of  estimating  clutter  statistics,  IFSTP  is  particularly 
suitable  for  applications  in  highly  non-homogeneous  or  unknown  clutter  environments. 

Target  detection  in  IFSTP  consists  of  two  steps.  The  first  step  is  2-D  Fourier 
transfonn  of  data  and  detecting  targets  and  interference  in  noise;  the  second  step  is 
extracting  targets  from  interference  based  on  their  different  features.  If  targets  and 
interference  are  clearly  separated  in  the  angle-Doppler  domain,  which  is  normally  true, 
separating  target  detection  from  interference  in  the  second  step  is  very  dependable  and 
the  overall  detection  probability  and  false-alann  rate  for  IFSTP  are  basically  determined 
by  processing  in  the  first  step,  i.e.,  target  detection  in  white  noise.  Based  on  the 
processing  described  in  the  first  step,  the  output  SNR  is  maximally  increased  NM  times 
the  input  SNR,  where  N  and  M  are  the  number  of  coherent  pulses  during  a  CPI  and  the 
number  of  antenna  elements,  respectively.  As  a  result,  the  detection  performance  of 
IFSTP  approaches  that  of  the  optimum  STAP  or  matched  filter,  provided  that  the  targets 
and  interference  are  clearly  separated  in  the  angle-Doppler  domain. 

Our  preliminary  analysis  also  indicates  that  effects  such  as  sensor  crabbing  and 
near-field  scattering,  which  might  affect  STAP  significantly,  will  not  affect  IFSTP 
performance.  Although  a  low  PRF  is  used  in  our  examples,  wavefonns  with  high  PRF  or 
medium  PRF  are  more  desirable  for  IFSTP  because  they  are  more  likely  to  separate 
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targets  further  from  interferes  in  the  angle-Doppler  plane.  However,  IFSTP  is  not  meant 
to  replace  traditional  STAP  or  any  other  adaptive  algorithm.  One  of  obvious 
shortcomings  of  IFSTP  is  that  the  numbers  of  array  elements  and  coherent  pulses  during 
a  CPI  need  to  be  large  enough  to  formulate  a  meaningful  image  in  the  angle-Doppler 
plane.  However,  this  problem  could  be  addressed  using  super-resolution  techniques  [9]. 
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PART  II:  Concurrent  Extraction  of  Target  Range  and 
Doppler  Information  by  Using  Orthogonal  Coding 
Waveforms 

Summary 

In  this  paper,  a  novel  approach  for  concurrent  extraction  of  target  range  and  Doppler 
information  in  radar  systems  is  introduced.  A  pair  of  orthogonal  waveforms  with 
different  pulse  repetition  frequencies  (PRF)  are  bundled  and  transmitted  as  a  single 
hybrid  radar  waveform  for  target  measurement.  At  the  receiver,  two  orthogonal  signal 
processing  channels  consisting  of  a  low  PRF  and  a  high  PRF  channel,  perform  the 
extraction  of  target  range  and  Doppler  information,  respectively.  The  CLEAN  algorithm 
is  innovatively  used  to  suppress  the  interference  from  the  cross-correlation  residues  and 
the  autocorrelation  sidelobe  in  the  LPRF  channel.  The  processing  results,  based  on  the 
proposed  signal  processing  approach,  demonstrate  the  technical  feasibility  of 
implementing  such  a  signal  processing  scheme  in  radar  systems. 48 
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1.  Introduction 

The  choice  of  the  radar  pulse  repetition  frequency  (PRF)  significantly  affects  the 
radar  performance,  relative  to  the  extraction  of  target  information  [1-3].  Traditionally, 
radar  systems  transmit  and  process  waveforms  with  fixed  PRFs,  which  may  be  high, 
medium  or  low,  during  a  coherent  processing  interval  (CPI).  If  a  low  PRF  (LPRF) 
waveform  is  used,  the  target  range  infonnation  is  instantly  acquired  though  time  delay 
measurement;  while  the  Doppler  frequency,  which  is  readily  converted  to  the  radial  target 
velocity  due  to  the  target  radial  motion,  is  normally  ambiguous  and  not  directly 
measurable.  Conversely,  if  a  high  PRF  (HPRF)  waveform  is  used,  the  target  Doppler 
shift  may  be  obtained  directly  through  a  Fourier  transform  of  the  echoes,  but  the  target 
range  becomes  ambiguous  and  has  to  be  indirectly  calculated  by  transmitting  multiple 
staggered  PRF  wavefonns  [4,  5].  Thus,  for  current  radar  systems,  if  a  waveform  with  a 
particular  PRF  is  used  for  direct  measurement  of  either  target  range  or  Doppler  shift,  the 
measurement  of  the  other  target  information  is  unavoidably  indirect,  and  is  thus  less 
accurate  and  more  time-consuming  [6,  7].  In  other  words,  existing  radar  systems  do  not 
have  the  capability  to  measure  target  range  and  velocity  simultaneously  and  accurately. 
However,  for  some  applications  such  as  tracking  high-speed  maneuvering  targets,  it  is 
highly  desirable  to  acquire  both  target  range  and  velocity  simultaneously  and  instantly  [8- 
10].  In  this  work,  we  will  demonstrate  the  technical  feasibility  of  implementing  such  a 
radar  system  by  transmitting  and  processing  hybrid  radar  waveforms  consisting  of 
orthogonal  coding  signals,  using  both  HPRF  and  LPRF  signals. 
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2.  Design  and  processing  of  hybrid  coding  waveform 

2.1.  Orthogonal  coding  waveforms 

The  hybrid  coding  waveform  consists  of  two  different  coded  sub-waveforms  sx  (?)  and 
s2  (?)  with  one  of  them  transmitted  in  HPRF  and  the  other  in  LPRF.  The  two  sub¬ 
waveforms  should  be  orthogonal  to  each  other  so  that  they  are  separable  through 
correlation  processing  at  the  receiver  for  target  range/Doppler  information  extraction,  i.e. 
they  should  satisfy  [11,  12]: 

jtfj  [t-t ^jdt » 0  VteR  (7) 

? 

For  the  transmission  and  processing  of  the  hybrid  wavefonn  to  be  compatible 
with  current  radar  system,  the  two  orthogonal  coding  wavefonns  are  interleaved,  bundled 
and  transmitted  as  a  single  radar  signal,  as  shown  in  Figure  8.  The  sub-waveform  sx  (?)  in 
the  figure  is  transmitted  with  low  PRF  for  instant  measurement  of  target  range.  The  sub¬ 
waveform  s2  (?)  is  transmitted  with  high  PRF  for  direct  measurement  of  target  Doppler 

frequency.  The  multiple  pulses  of  the  coding  waveform  s2  (?)  are  fully  phase-coherent 
for  pulse  Doppler  processing  at  the  receiver.  The  difference  between  the  coding 
waveform  s2  (?)  and  the  coherent  pulse  train  used  in  regular  pulse  Doppler  radar  is  that 

each  sub-pulse  in  s2  (?)  is  further  phase-coded  with  a  code  length  of  N. 
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Figure  8:  Orthogonal  coding  waveform  containing  both  HPRF  and  LPRF  sub-waveforms 


Assuming  that  the  coding  phase  sequences  for  the  coding  waveforms  s,  (Y)  and 
s2(t )  are  (l),$  (2),K  (A)]  and  j^2  (l),^2  (2) ,K  ,^2  ( A)}  ,  respectively,  where  A 

is  the  code  length,  the  coding  waveforms  Sj  (Y)  and  s2  (Y)  can  be  expressed  as 


h  (0  =  Z  R  (*  ~  nT )  exP  [M  («)] 


(8) 


and 


!  (0  =  Z  R  {*  - m )  eXP  [m  («)]  > 


(9) 


where  r  is  the  subpulse  width  of  the  coding  waveforms  and  R(t)  is  a  rectangular  pulse 


defined  as: 


*(<)= 


[1  It  0<t<r 

I  0  t  >  z 


(10) 


Both  wavefonns  sl  (Y)  and  s2  (Y)  are  assumed  to  be  of  pulse  duration  T,  where 


T  =  Ar . 


(ID 
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The  autocorrelation  functions  of  the  waveforms  sx  (/  )  and  s2  (/  )  in  (2)  and  (3)  and  their 
cross-correlation  function  are  assumed  to  be  Ax{t),  A2(t),  and  C(7),  respectively. 
Because  the  two  coding  signals  are  orthogonal  and  possess  pseudo-noise  characteristic, 


Ax  (7) ,  A2  (7) ,  and  C  (7)  may  be  approximated  as: 

(12) 

A2(t)=\s2(t')s*2(t'+t)dt' , 

(13) 

C(7)  =  J (t')s*  {t'+t^dt' xQ . 

(14) 

t' 


The  orthogonal  wavefonn  containing  sub-waveforms  y  (7)  and  ,s-2  ( /  ) ,  as  shown  in  Fig. 
8,  is  expressed  as  follows: 

H(t)=sx(t)  +  j^s2(t~T-l.T2)  (15) 

/= 1 

where  L  is  the  number  of  coherent  pulses  during  a  CPI  for  the  waveform  s2  (l )  and  7?  is 
the  pulse  repetition  period  for  the  s2  (/ )  waveform.  Because  y  (/  )  is  orthogonal  to 
s2  (7) ,  the  target  information  conveyed  in  the  two  signals  can  be  easily  separated  and 
extracted  at  the  receiver  through  a  correlation  processing. 
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2.2.  Processing  of  Orthogonal  Waveforms  for  Concurrent  Target 
Range  and  Doppler  Extraction 

If  a  received  radar  signal  contains  the  echoes  from  K  targets  during  a  CPI  at  an 
antenna  beam  dwelling,  the  amplitudes,  initial  phases,  time  delays  and  Doppler  shifts 
area  assumed  to  be  Mk , </>k ,  rkd  ,  and  fd  ,k  =  1, 2,L  , K ,  respectively.  The  received  signal 
prior  to  the  processing  is: 


£(0=2>, 

k= 1 


ej<t,k  s 


i  (f-Td) 


eJ2  «f!‘+Y^MkeJ 

k= 1  /= 1 


t-T-r l- 


I.T2  )e 


(16) 


Because  the  Doppler  shift  is  normally  much  smaller  than  the  reciprocal,  i.e.  \/T,  of  the 
duration  of  the  coding  signals  Sj(7)  or  s2(t),  its  effect  on  the  correlation  processing  in 
Fig.  9  is  ignored. 
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Figure  9:  Signal  processing  scheme  for  concurrent  extraction  of  target  ranges  and 

Doppler  frequencies 


The  target  range  infonnation,  conveyed  by  the  low  PRF  signal  .st,  (/ )  is  acquired  by 


correlating  the  signal  in  (16)  with  a  correlator  matched  to  s{  (7) .  This  correlator  is  termed 


as  the  low-PRF  (LPRF)  processor,  and  its  output  is  given  by: 


21 


PLr ( t ) * A, r‘)elw1'  +  Xi«,e*C(/-r- v‘ -I ■  I \Y”'“  ■ 


(17) 


k=\ 


k= 1  /= 1 


Because  sx  (*)  and  s2  ( / )  are  orthogonal  to  each  other,  the  magnitude  of  the  output  is 
approximately  equal  to: 


k= 1 


(18) 


The  ranges  of  the  radar  targets  are  calculated  from  the  delays  of  the  pulses  in  (12)  from 
the  moment  when  the  wavefonn  is  transmitted. 


Similarly,  a  correlator  matched  to  5,  (/ )  is  employed  in  the  high  PRF  channel  to  extract 
target  Doppler  information  from  the  received  signal  in  (16).  The  correlation  result  is: 

P„n (0 - Z Mke^c(t - r) )e-/2^'  +  X Z AVM A2{t -T -zkd -l ■  T2 )e>2^‘  •  (19) 

k= 1  k= 1  /= 1 


Because  the  magnitude  of  the  cross-correlation  function  C  is  small  and  the 
autocorrelation  function  Ai  closely  approximates  the  Dirac  function,  the  above  equation 
may  be  further  simplified  as: 


nxfiiT2 


(20) 


where  0;  =(j)k+2 n //  ( T  +  xkd  j  is  a  fixed  value 


Generally,  the  output  signal  in  (20)  contains  multiple  targets  locating  at  different 
range  cells  (i.e.,  time  delays).  To  acquire  Doppler  information  of  the  k- th  target  for  k=  1, 
2,  ...,  K,  the  signal  in  (20)  is  sampled  at  t  =  T  +  rkd+l-T2  ;  /  =  1,2,K  ,L,  which  can  be 
detennined  using  the  target  range  obtained  from  (12).  The  Fourier  transform  is 
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performed  on  the  L  sampled  data  from  (20)  for  each  target.  If  the  received  signal  energy 
loss  due  to  transmitting  eclipse  is  ignored,  the  transform  result  is  approximately  equal  to: 

|F<(/)|=!I-Mi-A(0)5(/-/;);-l/(2r2)</<l/(2r;);t  =  l,2,L  ,K  (21) 

where  8  is  the  Dirac  function  in  the  frequency  domain.  Noticeably,  the  processing  result 
for  Doppler  information  in  (17)  is  much  stronger  than  that  for  the  range  information 
shown  in  (14)  due  to  the  two-time  coherent  integrations,  i.e.,  correlation  and  Fourier 
transfonn. 

For  the  two-channel  processing  results  of  multiple  target  signals  in  (18)  and  (21),  the 
range  and  Doppler  data  of  the  same  target  are  separately  obtained,  and  thus  a  pairing 
processing  of  range  and  Doppler  infonnation  is  needed  for  each  of  the  targets.  For  the 
high  PRF  channel  processing,  the  range  information  in  (20)  is  ambiguous.  A  simple 
approach  is  to  use  the  actual  range  of  a  target  in  (18)  to  calculate  its  ambiguous  range  and 
perform  a  Fourier  transform  only  on  the  data  sampled  from  the  ambiguous  range  unit.  An 
alternative  method  is  to  perform  the  Fourier  transfonn  on  all  ambiguous  range  units  in 
(20)  for  the  high-PRF  channel  and  to  detect  all  targets  with  their  Doppler  frequencies  and 
ambiguous  ranges.  If  a  target’s  Doppler  and  ambiguous  range  are  f D  and  R , , 
respectively,  its  actual  range  R  measured  using  the  low  PRF  channel,  is  associated  with 
the  Doppler  frequency  fD  by  the  following  equation: 


where  C  is  the  speed  of  light  and  _x\  denotes  the  largest  integer  not  greater  than  x. 
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Apparently  from  (18)  and  (21),  the  ranges  and  the  Doppler  frequencies  of  multiple  targets 
are  directly  and  simultaneously  measurable  by  transmitting  and  processing  the  orthogonal 
waveform  scheme  proposed  in  this  paper.  The  signal  processing  scheme  for  the 
concurrent  range  and  Doppler  measurement  is  shown  in  Fig.  9. 

2.3.  Interference  reduction  processing  using  the  CLEAN 
algorithm 

The  results  in  (17)  and  (19)  out  of  the  correlation  processing  in  the  HPRF  and  LPRF 
channels  the  cross-correlation  function  C(  t  )  of  two  orthogonal  coding  signals  y  (/ )  and 

s2  (t)  is  assumed  to  be  zero  and  ignored.  Because  the  actual  cross-correlation  function  is 

not  zero,  in  the  multi-target  detection  case,  the  cross-correlation  residues  from  strong 
targets  could  be  comparable  to  the  responses  of  weak  targets,  resulting  in  detection 
degradation  of  the  weak  targets.  For  the  same  reason,  the  sidelobe  of  a  strong  target 
response  could  interfere  with  the  detection  of  a  weak  target.  Therefore,  the  CLEAN 
algorithm  is  used  to  reduce  the  cross-correlation  resides  between  the  HPRF  and  LPRF 
waveforms  as  well  as  the  sidelobes  of  strong  target  responses.  However,  for  the  HPRF 
channel  processing,  the  target  signal  is  coherently  integrated  twice;  as  a  result,  even  the 
response  of  a  weak  target  is  normally  much  stronger  than  the  sidelobe  or  cross¬ 
correlation  residues  from  a  strong  target.  Therefore,  interference  reduction  is  only 
considered  for  the  range  infonnation  extraction  processing  (i.e.,  Low  PRL  channel). 

The  CLEAN  algorithm  has  been  widely  used  to  remove  interference  from  the 
autocorrelation  sidelobe,  using  matched  processing  in  the  case  of  weak  target  detection  or 
clear  image  formulation  [13-17].  In  this  work,  the  CLEAN  algorithm  is  adapted  to 
remove  both  the  autocorrelation  sidelobe  and  the  cross-correlation  residues  in  the  low- 
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PRF  channel.  The  critical  part  of  implementing  the  algorithm  is  to  locate  the  positions  of 
the  cross-correlation  functions  and  to  determine  their  magnitudes  and  phases.  The  process 
is  to  start  from  the  correlation  result  in  the  low-PRF  channel.  For  each  of  the  detected 
targets,  there  are  L  (number  of  the  coded  pulses  for  the  high-PRF  channel)  cross¬ 
correlation  functions  that  need  to  be  deducted  from  the  result.  The  location  of  the  cross¬ 
correlation  functions  can  be  determined  from  the  detected  target  location  in  the  low-PRF 
channel  and  the  relative  time  delays  from  the  low-RPF  pulse  to  the  high-PRF  pulses.  The 
magnitude  can  be  directly  derived  from  the  magnitude  of  the  mainlobe  of  the  target 
response,  while  the  phase  needs  to  be  estimated  from  the  target  mainlobe  phase  and  the 
target  Doppler  frequency  measured  in  the  high-PRF  channel.  The  procedure  for 
implementing  the  proposed  interference  reduction  method  is  described  below: 

STEP  1:  If  the  low-PRF  correction  processing  result  is  a  data  sequence 
jc/(z)  /  =  1, 2, 3,K  | ,  the  strongest  target  is  detected  as  the  index  J  with  magnitude  M  j , 

phase  y/j ,  and  Doppler  frequency  fJD  (measured  at  the  high-PRF  channel),  the  target 

sidelobe  sequence  |^4  (n);n=l,2,K , 2Af — 1|  is  deducted  from  the  data 

sequence  {d(i)}  with  A  (AT)  aligned  to  d (J ) ,  where  the  data  sequence  is  derived  from  the 
autocorrelation  function  of  the  low-PRF  coding  signal: 


A'(n)  = 


M  pJ¥j 

J  A ,  (/;)  ;  N 


N 


;  n  =  N 


(23) 
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STEP  2:  Subtract  the  cross-correlation  sequence  jc,  (n)  ;  n  =  l, 2, K  ,27V-l|  between 
the  low-PRF  coding  pulse  and  the  /-th  high-PRF  coding  pulse  (1  <  /  <  L)  from  the  data 
sequence  {d  (/)}  for  the  target.  The  cross-correlation  sequence  is: 


C\  (n) 


M  eAvj+2^fDj-^i ) 

N 


C(n) 


(24) 


where  Ar,  is  the  time  delay  from  the  low-PRF  coding  pulse  to  the  /-th  high-PRF  coding 
pulse  of  the  orthogonal  signal,  and  C(n)  is  the  cross-correction  function  of  s^t)  and 
s2(f). 

STEP  3:  Repeat  the  processing  in  STEP  2  until  all  cross-correlation  resides  are  removed 
from  the  data  sequence,  i.e.  /  =  1,2,L  ,L  . 

STEP  4:  Find  the  next  largest  target  signal  from  the  data  sequence  and  repeat  Steps  1-3 
until  no  detectable  target  signal  exists  in  the  data  sequence.  The  final  data 
sequence  {d(i);i  =1,2,3, K}  is  absent  of  the  autocorrelation  sidelobe  or  the  cross¬ 
correlation  residues. 

The  uniqueness  of  the  CLEAN  algorithm  as  applied  in  this  work  is  to  subtract  the 
interference  from  multiple  cross-correlation  functions  rather  than  the  sidelobes  of  the 
autocorrelation  functions.  In  addition,  to  correctly  implement  the  CLEAN  algorithm,  one 
needs  to  estimate  and  use  the  target  Doppler  frequency  and  the  initial  phase  of  the 
target’s  echo.  Like  any  other  CLEAN  algorithm,  the  processing  result  is  slightly  tainted 
by  the  inaccuracy  of  the  magnitudes  of  strong  target  signals  that  are  affected  by  additive 
white  noise  and  the  overlapping  sidelobe  of  other  nearby  targets  [17].  Since  the 
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magnitude  of  the  strong  target  signals  is  normally  much  larger  than  that  of  white  noise  or 
other  target  sidelobe,  the  target  magnitude  distortion  is  marginal.  In  this  application,  the 
CLEAN  result  is  also  inconsequentially  affected  by  the  accuracy  of  the  target  Doppler 
frequency  estimation  in  the  high-PRF  channel. 

3.  Processing  Results 

The  proposed  radar  signal  processing  is  applied  to  the  radar  data  simulated  from  a 
scenario  with  multiple  moving  targets.  The  orthogonal  waveform  is  designed  using 
polyphase  coding  sequences  by  simulated  annealing  (SA)  [10].  The  code  length  N  is  64 

and  the  admissible  phase  values  for  the  code  design  are  =  0,1,2,K  ,7j  .  The 

relevant  parameters  of  the  designed  orthogonal  waveform  shown  in  Fig. 8  are: 

•  Pulse  period  of  the  LPRF  waveform  s,  (t) :  T\  =  8.0ms 

•  Pulse  period  of  HPRF  waveform  s2  (t) :  T2  =  90ps 

•  Pulse  duration:  T=  lOps 

•  Subpulse  duration:  r  =  TIN  =  (5/32)ps 

•  Number  of  pulses  during  a  CPI  for  the  HPRF  waveform  s2  (t)  :  L  =  80 

It  is  assumed  that  the  Doppler  frequency  fD  and  the  coding  signal  pulse  duration  meet 
the  following  requirement  [18]: 

fo-T«  1  (25) 
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3.1.  Processing  results  without  interference  reduction 

If  the  received  radar  data  contains  a  single  moving  target  at  a  range  of  150km  and 
Doppler  frequency  of  -2.5kHz,  the  processing  results  of  the  LPRF  and  HPRF  channels 
are  shown  in  Figs  12  and  11,  respectively.  The  time  delay  of  the  target  pulse  in  Fig.  10 
from  the  transmitting  time,  i.e.  at  time  zero,  exactly  reflects  the  target  range.  Besides  the 
target  signal,  the  interference  from  the  cross-correlation  between  the  LPRF  channel  signal 
Sj  (t)  and  the  HPRF  channel  signal  s2  ( / )  is  quite  visible  in  the  result.  This  kind  of 

interference  signal  might  mask  weak  targets  and  thus  preventing  them  from  being 
detected  in  the  multi-target  case. 


Figure  10:  Normalized  output  of  the  LPRF  processing  channel  with  a  single  target  at  a 

range  of  150km 
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For  the  FIPRF  channel  processing,  Fig.  11  shows  the  result  of  Doppler  extraction 
only,  based  on  the  data  sampled  at  the  target  range  cell.  The  Fourier  transform  results  at 
other  range  cells  are  below  the  target  detection  threshold.  The  frequency  in  Fig.  1 1  is 
exactly  the  target’s  Doppler  frequency  -2.5Hz.  Due  to  the  two-time  coherent  integrations 
for  target  signal  in  the  FIPRF  channel,  as  indicated  in  (21),  the  mainlobe  of  the  target 
signal  is  much  stronger  than  the  sidelobe  of  the  target  response  and  the  cross-correlation 
residues.  As  a  result,  the  interference  is  almost  indiscernible  for  the  FIPRF  channel 
processing  result. 


-2  0  2 

Doppler  Shift  (kHz) 

Figure  11:  Normalized  output  of  the  HPRF  processing  channel  for  the  same  radar  data  in 
Fig.  10  containing  a  single  target  with  Doppler  frequency  of  -2.5kHz. 
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For  practical  radar  operation,  it  is  most  likely  that  multiple  targets  exist  in  the 


same  antenna  beam  position.  Flence,  a  simulated  radar  echo  containing  three  moving 


targets  is  applied  to  the  proposed  signal  processing  scheme.  The  ranges  of  the  three 


targets  are  60km,  170km  and  350km;  their  relative  magnitudes  are  1,  0.8,  and  0.5;  and 


their  radial  Doppler  frequencies  are  1kHz,  3kHz  and  -1kHz,  respectively.  The  additive 


white  noise  with  SNR  equal  to  14dB  is  added  to  the  radar  data.  The  low-PRF  channel 


processing  result  for  the  multi-target  radar  signal  is  shown  in  Fig.  12.  The  delayed  times 


for  the  three  targets  shown  in  the  figure  correspond  to  their  actual  ranges.  The  cross¬ 


correlation  residues  from  the  HPRF  waveform  are  also  observable  in  the  result. 
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Figure  12:  LPRF  channel  processing  result  for  three  targets  at  ranges  60km,  170km,  and 
350km  and  magnitudes  of  1,  0.8,  and  0.5,  respectively.  Additive  white  noise  with  SNR=14  is 

included  in  the  input  data. 
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Theoretically,  for  the  HPRF  channel  processing,  the  Fourier  transfonn  or  Fast 
Fourier  Transform  (FFT)  needs  to  be  perfonned  at  all  ambiguous  range  cells  during  a 
FIPRF  CPI,  and  subsequently  target  detection  is  performed  at  all  ambiguous  ranges.  The 
processing  results,  i.e.,  the  target  Doppler  frequencies,  for  the  ambiguous  range  cells  in 
which  the  maximum  output  value  is  above  the  detection  threshold  are  shown  in  Fig  13. 


Figure  13:  HPRF  processing  result  for  the  same  radar  data  in  Fig.  12  with  target  Doppler 

frequencies  of  1kHz,  3kHz  and  -1kHz. 

The  white  noise  in  the  input  data  has  slightly  increased  the  interference  level  of 
the  processing  results  in  the  HPRF  and  LPRF  channels,  as  shown  in  Figs  12  and  13,  but 
the  overall  system  performance  is  not  impacted  by  the  additive  white  noise  or  any  other 
inference  signal  in  the  input  data. 
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The  remaining  issue  in  the  multiple  target  detection  case  is  how  to  associate  a 
range  measured  in  the  LPRF  channel  with  a  Doppler  frequency  measured  in  the  HPRF 
channel  for  the  same  target.  A  straightforward  approach  is  to  measure  the  ranges  and  the 
Doppler  frequencies  separately  and  independently;  and  to  pair  the  range  and  Doppler  for 
each  of  the  targets  using  (22).  The  alternative  approach  is  to  use  the  measured  range  for  a 
target  from  the  LPRF  channel  to  calculate  its  ambiguous  range  and  perform  Doppler 
processing  only  at  the  ambiguous  range  cells.  The  latter  approach  is  simple,  but  slightly 
delays  the  execution  of  Doppler  information  extraction.  Using  either  approach,  one  can 
conveniently  find  that  the  range-Doppler  pairs  from  Figs  12  and  13  for  the  three  targets  to 
be:  (60km,  1kHz),  (170km,  3kHz)  and  (350km,  -1kHz).  It  is  possible  that  two  or  more 
targets  may  actually  have  the  same  ambiguous  range,  which  effectively  invalidates  either 
of  the  two  described  approaches  for  pairing  the  range  and  the  Doppler  frequency  of  a 
target  among  multiple  targets.  In  this  case,  a  slightly  different  PRF  has  to  be  used  for  the 
HPRF  channel  to  resolve  the  ambiguous  ranges  of  two  or  more  targets  with  different 
ranges. 

3.2.  Processing  result  with  the  CLEAN  algorithm  applied 

In  the  multiple  target  scenarios,  as  discussed  in  Section  II,  the  detection  of  the 
weak  target  signals  in  the  low-PRF  channel  could  be  interfered  by  the  sidelobe  of  the 
stronger  target  responses  and  the  cross-correlation  residues  between  the  low-PRF  channel 
correlation  filter  and  the  high-PRF  coding  signals  reflected  by  stronger  targets.  Fig.  14 
shows  the  LPRF  channel  processing  result  for  a  radar  signal  containing  targets  with 
relative  signal  magnitudes  of  1,  0.8,  0.2  and  0.12.  The  target  ranges  are  50km,  120km, 
230km  and  300km;  and  their  respective  Doppler  frequencies  are  4kHz,  -3kHz,  2kHz,  and 
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0.5kHz.  With  both  strong  and  weak  targets  existing  in  the  same  radar  echo,  the  two  weak 
target  signals  are  completely  masked  by  the  sidelobe  of  the  strong  target  signals  and  the 
cross-correlation  residues  between  the  LPRF  filter  and  the  HPRF  waveforms.  Therefore, 
the  interference  signal  needs  to  be  completely  removed  or  at  least  significantly  reduced  to 
effectively  detect  weak  targets.  In  this  work,  the  CLEAN  algorithm  is  applied  to  reduce 
the  interference  from  the  LPRF  channel  processing  result.  With  the  interference- 
suppressed  processing,  the  ranges  of  the  four  targets  including  the  weak  ones  are 
correctly  displayed  in  Fig.  15.  The  target  Doppler  frequencies,  which  are  shown  in  Fig. 
16,  are  required  for  correctly  implementing  the  CLEAN  algorithm.  Note  that  the 
CLEAN  algorithm  does  not  remove  the  interference  completely  because  the  correlation 
processing,  prior  to  the  CLEAN  processing,  has  slightly  distorted  the  pulse-to-pulse 
phase  shifts  caused  by  the  Doppler  frequency.  In  addition,  existence  of  white  noise  in  the 
received  signal  also  makes  perfect  interference  removal  impossible. 
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Figure  14:  LPRF  detection  result  for  four  targets  with  magnitudes  of  1,  0.8,  0.2,  and  0.12. 
The  two  weakest  targets  cannot  be  detected  due  to  the  interference  from  cross-correlation 

residues. 
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Figure  15:  LPRF  channel  processing  result  for  the  same  radar  data  in  Fig.  14.  The  target 
ranges  are  50km,  120km,  230km  and  300km.  The  two  weak  targets  become  detectable  after 
interference  suppression  using  the  CLEAN  algorithm. 
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Figure  16:  Measured  Doppler  frequencies  of  the  four  targets  for  the  data  used  in  Fig.  14. 
The  extracted  Doppler  frequencies  are  4kHz,  -3kHz,  2kHz,  and  0.5kHz.  The  Doppler 
extraction  is  performed  after  interference  reduction  processing 


Figure  17  displays  the  interference  reduction  processing  for  the  same  input  data 
used  in  the  processing  of  Fig.  12,  which  is  corrupted  by  additive  white  noise.  The 
processing  result  indicates  that  the  effectiveness  of  the  interference  reduction  algorithm  is 
not  affected  by  the  additive  unknown  white  noise. 
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Figure  17:  The  CLEAN  algorithm  is  applied  to  the  LPRF  channel  processing  result  for  the 
same  input  data  as  in  Fig.  13  with  additive  white  noise  in  the  input  data 


The  computational  complexity  of  implementing  the  CLEAN  algorithm  in  this 
work  is  about  Q  L-  2  •  (2 N  - 1)  additions  and  multiplications,  where  Q  is  the  number  of 
targets  in  a  CPI.  Because  L  and  N  are  fixed  values  and  Q  is  not  normally  very  large  for 
point-target  model,  the  computational  cost  for  using  the  CLEAN  algorithm  is  not 
expensive. 
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4.  Conclusions 

Theoretical  analysis  and  simulation  results  have  demonstrated  that  the  ranges  and 
velocities  of  high-speed  targets  can  be  concurrently  measured  by  transmitting  and 
processing  a  specially  designed  orthogonal  waveform  scheme.  Results  show  that  the 
proposed  signal  processing  system  performs  the  desired  function  of  concurrent 
measurement  exactly  as  expected.  The  application  of  the  CLEAN  algorithm  is  necessary 
to  suppress  interference  in  the  LPRF  channel  for  reliable  detection  of  weak  targets. 

The  proposed  orthogonal  waveform  processing  scheme  is  compatible  with 
traditional  radar  signal  processing  techniques  for  clutter  rejection  such  as  Moving  Target 
Indication  in  LPRF  channel  [2]  or  Pulse  Doppler  processing  in  HPRF  channel  [7].  They 
may  be  directly  cascaded  to  or  combined  into  the  proposed  signal  processing  scheme 
without  degrading  performance.  Therefore,  the  existence  of  clutter  in  the  received  radar 
signal  does  not  impair  the  generality  of  the  proposed  algorithm. 

The  new  technique  may  be  applied  for  effective  detection  and  tracking  of  high¬ 
speed  maneuvering  targets  such  as  high-speed  aircraft  and  ballistic  missiles  because  of 
the  availability  of  target  range  and  radial  velocity  once  a  target  is  detected.  Similarly, 
more  concurrent  radar  functions  such  as  target  imaging  or  high-resolution  range  profiling 
may  be  further  integrated  by  including  additional  orthogonal  waveform  such  as  a 
wideband  waveform  as  part  of  the  transmitting  signal  and  the  corresponding  processing 
channel  at  the  receiver.  The  newly  added  waveform  must  be  orthogonal  to  both  existing 
coding  waveforms. 

For  the  HPRF  channel  in  the  proposed  signal  processing  scheme,  due  to  the  two- 
time  integrations,  the  SNR  of  the  output  signal  in  this  channel  is  much  higher  than  that  in 
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the  LPRF  channel.  If  the  coding  length  of  the  LPRF  waveform  is  designed  to  be  larger, 
the  SNRs  of  the  two  channels  will  be  more  balanced.  Therefore,  it  is  more  desirable  to 
design  a  pair  of  orthogonal  waveforms  with  different  coding  lengths  in  the  future. 
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